All computer codes are in MATLAB. A brief description of the codes is as follows:

-'Bloch_phonon.m': Bloch wave calculation of 'zero' loss energy filtered diffraction intensities. Includes only phonon inelastic scattering. 
-'Bloch_single_plasmon_phonon.m': Bloch wave calculation of single plasmon energy filtered diffraction intensities. Includes multiple phonon scattering and single plasmon scattering. 
-'Bloch_double_plasmon_phonon.m': Bloch wave calculation of double plasmon energy filtered diffraction intensities. Includes multiple phonon scattering and double plasmon scattering. 

-'Bloch_phonon_DP.m': Bloch wave calculation of thermal diffuse scattering intensity distribution (diffraction pattern, DP) for single phonon scattering.
-'Bloch_plasmon_DP.m': Bloch wave calculation of plasmon diffuse scattering intensity distribution (diffraction pattern, DP) for single plasmon scattering (no phonons).

All of the above codes are for [110]-Si.